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ABSTRACT 

J> ! We perform a systematic study of the one-loop renormalization-group-improved effective 
^ . potential of the minimal supersymmetric extension of the Standard Model (MSSM), includ- 
ing CP violation induced radiatively by soft trilinear interactions related to squarks of the 
third generation. We calculate the charged and neutral Higgs-boson masses and couplings, 
including the two-loop logarithmic corrections that arise from QCD effects, as well as those 
associated with the top- and bottom-quark Yukawa couplings. We also include the poten- 
tially large two-loop non-logarithmic corrections induced by one-loop threshold effects on 
the top- and bottom-quark Yukawa couplings, due to the decoupling of the third-generation 
squarks. Within this minimal CP-violating framework, the charged and neutral Higgs sec- 
tors become intimately related to one another and therefore require a unified treatment. 
In the limit of a large charged Higgs-boson mass, M H ± ^> Mz, the lightest neutral Higgs 
boson resembles that in the Standard Model (SM), and CP violation occurs only in the 
heavy Higgs sector. Our analysis shows that sizeable radiative effects of CP violation in 
the Higgs sector of the MSSM may lead to significant modifications of previous studies for 
Higgs-boson searches at LEP 2, the Tevatron and the LHC In particular, CP violation 
could enable a relatively light Higgs boson to escape detection at LEP 2. 
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1 Introduction 



The violation of CP was first observed in the neutral kaon system [|], where a deviation 
from the superweak theory has recently been confirmed ||, and CP violation in B- meson 
decays || is strongly suggested by recent experiments ||]. In addition to its interest for 
particle physics, CP non- conservation provides a key ingredient for cosmological baryoge- 
nesis, namely for explaining the underlying mechanism which caused matter to dominate 
over anti-matter in our observable universe |J . Although a fundamental understanding of 
the origin of CP violation is still lacking, most of the scenarios proposed in the existing 
literature indicate that Higgs interactions play a key role in mediating CP violation. For 
instance, CP violation is broken explicitly in the Standard Model (SM) by complex Yukawa 
couplings of the Higgs boson to quarks. Another appealing scheme of CP violation occurs 
in models with an extended Higgs sector, in which the CP symmetry of the theory is bro- 
ken spontaneously by the ground state of the Higgs potential ||. Supersymmetric (SUSY) 
theories, including the minimal supersymmetric extension of the Standard Model (MSSM), 
predict an extended Higgs sector, and therefore may realize either or both the above two 
schemes of explicit and spontaneous CP violation. However, the Higgs potential of the 
MSSM is invariant under CP at the tree level, and any explicit or spontaneous breakdown 
of CP symmetry can arise only via radiative corrections. The case of purely spontaneous 
CP violation in the MSSM J7J leads to an unacceptably light CP-odd scalar ||, as a result 
of the Georgi-Pais theorem |J, and hence such a scenario is ruled out experimentally. 

It has recently been shown that the tree-level CP invariance of the MSSM Higgs 
potential may be violated sizeably by loop effects involving soft CP-violating trilinear in- 
teractions of the Higgs bosons to top and bottom squarks. A detailed study |T| has shown 
that significant CP-violating effects of level crossing in the Higgs sector can take place in 
such a minimal SUSY scenario of explicit CP violation, which may lead to drastic mod- 



ifications of the tree-level couplings of the Higgs particles to fermions [|lT||12l and to the 
W and Z bosons fT^j. The latter can have important phenomenological consequences 
on the production rates of the lightest Higgs particle, even though the upper bound on 
its mass was found [[□]] to be very similar to that found previously in the CP-conserving 

case 



IJ[iJ,[IJ[lJ,|Tj,|T^.|2g.|2Tl. The MSSM predicts an upper bound on the lightest Higgs 



boson mass of approximately 110 (130) GeV for small (large) values of the ratio of Higgs 
vacuum expectation values tan/3 ~ 2 (> 15). On the other hand, experiments at LEP 2, 
running at center-of-mass energies yfs = 196-202 GeV, have placed a severe lower bound 
of approximately 108 GeV on the mass of the SM Higgs boson P2| . LEP 2 is expected to 
run at center-of-mass energies up to y/s = 206 GeV during the year 2000. Consequently, 
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for stop masses smaller than 1 TeV, a significant portion of the parameter space spanned 
by tan/3 and the CP-odd scalar mass can be tested for the CP-conserving case in this 



next round of experiments at LEP 2 [IT]. However, the explorable region of parameters is 
smaller for larger amounts of stop mixing and/or larger CP- violating phases. A decisive 
test of such a scenario can only be provided by the upgraded Tevatron collider and the 
LHC. 

The earlier study of the renormalization-group (RG) improved effective potential of 
the MSSM with explicit CP violation in ]TTJ] was based on an expansion of the Higgs quartic 
couplings in inverse powers of the arithmetic average of stop and sbottom masses, under 
the assumption that the mass splittings of the left- and right-handed stop and sbottom 
masses are small. The mass expansion of the one-loop effective potential was truncated up 
to renormalizable operators of dimension four. Although the above approach captures the 
basic qualitative features of the underlying dynamics under study, it is known [T^JTTJ] that 
such a mass expansion has limitations when the third-generation squark mixing is large. 
Since the dominant CP-violating loop contributions to the effective Higgs-boson masses and 
mixing angles occur for large values of the third-generation squark-mixing parameters, it is 
necessary to provide a more complete one-loop computation of the effective MSSM Higgs 
potential with explicit radiative breaking of CP invariance, including non-renormalizable 
operators below the heavy scalar-quark scale and without resorting to any other kinematic 
approximations.^ 

To this end, we consider here the two-loop leading logarithms induced by top- and 
bottom-quark Yukawa couplings as well as those associated with QCD corrections, by means 
of RG methods. In the calculation of the RG-improved effective potential, we also include 
the leading logarithms generated by one-loop gaugino and higgsino quantum effects |14 



Finally, we implement the potentially large two-loop corrections that are induced by the 
one-loop stop and sbottom thresholds in the top- and bottom-quark Yukawa couplings, 
which may become particularly relevant in the large-tan (3 regime. On the basis of the RG- 
improved effective potential, we present predictions for the Higgs-boson mass spectrum and 
the effective Higgs couplings to fermions and gauge bosons. The results of the analysis are 
compared with those obtained in the CP-conserving case ]l8[|l^] and also with the results 
obtained by truncating the one-loop RG-improved effective potential up to renormalizable 
operators |TT| . 

In this analysis, it is important to consider the constraints on the low- 



*We recall that the diagrammatic computation of scalar-pseudoscalar transitions is the focus of |T3], 
whilst sbottom contributions and relevant D-term effects on the effective potential are not considered 
in HI. 
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energy CP-violating parameters of the MSSM that originate from experimental up- 



per limits on the electron [23| and neutron [24| electric dipole moments (EDMs) 
2B|,g3,g|,gD|,gD|,^T],^|,H,gl|,^|. Most of the EDM constraints affect the CP- violating 



couplings of the first two generations [^5[^]. Thus, making the first two generations of 
squarks rather heavy, much above the TeV scale [p7 , p8 |, is a possibility that can drastically 
reduce one-loop contributions to the neutron EDM, without suppressing the CP-violating 
phases of the theory. Another interesting possibility for avoiding any possible CP crisis 
in the MSSM is to arrange for cancellations among the different EDM terms, either at 
the level of short- distance diagrams or (for the neutron EDM) at the level of the 



strong-interaction matrix elements for operators with s, u and d quark flavours p0|j31 



Alternatively, one might require a specific form of non-universality in the soft trilinear 
Yukawa couplings P2],p3 • 



However, third-generation squarks can also give rise by themselves to observable ef- 



fects on the electron and neutron EDMs through the three-gluon operator ||33|| , through 
the effective coupling of the 'CP-odd' Higgs boson to the gauge bosons and through 
two- loop gaugino/higgsino- mediated EDM graphs P5| . These different EDM contributions 
of the third generation can also have different signs and add destructively to the electron 
and neutron EDMs. In our phenomenological discussion, we take into account the relevant 
EDM constraints related to the CP-violating parameters of the stop and sbottom sectors. 

This paper is organized as follows: in Section 2 we calculate the complete one-loop 
CP- violating effective potential, and derive the analytic expressions for the effective charged 
and neutral Higgs-boson mass matrices. Technical details are given in Appendices A and 
B. Section 3 describes our approach to determining the RG-improved Higgs-boson mass 
matrices, after including the leading two-loop logarithms associated with Yukawa and QCD 
corrections. Section 4 is devoted to the calculation of the effective top- and bottom-quark 
Yukawa couplings, in which one-loop threshold effects of the third-generation squarks are 
implemented. In Section 5, we discuss the phenomenological implications of representative 
CP- violating scenarios compatible with EDM constraints for direct Higgs searches at LEP 2 
and the upgraded Tevatron collider. We also compare the results of our analysis with those 
obtained using the mass-expansion method [II]. Finally, in Section 6 we summarize our 
conclusions. 



2 CP- Violating One-Loop Effective Potential 

In this Section, we first describe the basic low-energy structure of the MSSM that contains 
explicit CP-violating sources, such as soft CP- violating trilinear interactions. Then we cal- 
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culate the general one-loop CP- violating effective potential. Finally, after implementing the 
minimization tadpole conditions related to the Higgs ground state, we derive the effective 
charged and neutral Higgs-boson mass matrices. 

CP violation is introduced into the MSSM through the Higgs superpotential and the 
soft supersymmetry-breaking Lagrangian: 



W = hiH{iT 2 LE + h d Hfir 2 QD + h u Q J ir 2 H 2 U - fiH(ir 2 H 2 , (2.1) 



- C 



soft 



-\(m~ g X1\1 + rriyy A-.A-. + mgAgAg + h.c.) + M\l)L + M%($Q 

+ M%U*U + M 2 D D*D + M'j iri-: + m?$t$i + m 2 2 ^ 2 - (Bfi^ir 2 ^ 2 
+ h.c.) + ( hiAi $\Ze + h d A d <&[QD - h u A u ^ir 2 QU + h.c), (2.2) 



where $1 = it 2 <&\ is the scalar component of the Higgs chiral superfield Hi and r 2 is 
the usual Pauli matrix. The conventions followed throughout this paper, including the 
quantum-number assignments of the fields under the SM gauge group, are displayed in 
Table § 



As can be seen from ( |2.1|) and (2.2), the MSSM includes additional complex pa- 



rameters with new CP-odd phases that are absent in the SM. These new CP-odd phases 
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Table 1: The field content of the MSSM. 
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may reside in the following parameters: (i) the mass parameter /i describing the bilinear 
mixing of the two Higgs chiral superfields in the superpotential; (ii) the soft supersymetry- 
breaking gaugino masses rrig, and of the gauge groups SU(3) C , SU(2)l and U(l)y, 
respectively; (iii) the soft bilinear Higgs-mixing mass Bfi; and (iv) the soft supersymmetry- 
breaking trilinear couplings Af of the Higgs bosons to sfermions. In addition, there may 
exist other large CP-odd phases, associated with flavor off-diagonal soft supersymmetry- 
breaking masses of squarks and sleptons. We assume that these off-diagonal masses are 
small and therefore do not give sizeable contributions to the effective Higgs potential |36[ . 



The number of independent CP-odd phases may be reduced if one prescribes a universality 
condition for all gaugino masses at the unification scale Mx', the gaugino masses will then 
have a common phase. Correspondingly, the different trilinear couplings Af may be con- 
sidered to be all equal at Mx, i.e., Af = A. In this case, however, because of the different 
RG running of the phases of the trilinear couplings, their values at low energies will be 
different |] Two CP-odd phases may further be eliminated by employing the following two 
global symmetries that govern the dimension-four operators in the MSSM Lagrangian: 

(i) The U(1) Q symmetry specified by Q{%) = 1, Q{%) = -2, Q(Q) = Q(L) = 0, 
Q(U) = 2 and Q(D) = Q(E) = —1. This U(1)q symmetry is broken by the \i 
parameter and the respective soft supersymmetry-breaking one, Bfi. 

(ii) The U(l)# symmetry acting on the Grassmann- valued coordinates, i.e., the 9 co- 
ordinate of superspace carries charge 1. Under the R transformation, the matter 
superfields and gaugino fields carry charge 1, whilst the Higgs superfields are R- 
neutral. The R symmetry is violated by the gaugino masses, the trilinear couplings 
Af and the parameter /i. 

We concentrate on the parameters which may have a dominant CP-violating effect 
on the MSSM Higgs potential, under the assumption of a common phase for the gauginos; 
the latter is made less important by the fact that the one-loop gaugino corrections are 
subdominant compared to the ones induced by the third-generation squarks. As has been 
mentioned above, two CP-odd phases of the complex parameters {/i, m\ 2 , m\, A} may 
be removed by employing the global symmetries (i) and (ii). Specifically, one of the Higgs 
doublets and the common phase of the gaugino fields can be rephased in a way such that 
the gaugino masses and Bfi become real numbers. As a consequence, arg(/x) and arg(A ttb ) 
are the only physical CP-violating phases in the MSSM which affect the Higgs sector in a 
relevant way. 



'''For a discussion of the RG effects, see 137 
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It is obvious from ( |2.1Q and ( |2.2| ) that the Yukawa interactions of the third-generation 
quarks, Q T = (t^, and t R , b R , as well as their SUSY bosonic counterparts, Q T = (t^, b L ) 
and t R , b R , play the most significant role in radiative corrections to the Higgs sector. 
Therefore, it is useful to give the interaction Lagrangians related to the F and D terms of 
the third generation: 



c 



D 



N 2 |$+Q| 2 + \h t \ 2 \^r 2 Q\ 2 
-(fih* b Q^ 2 b R + ^h* t Q ] iT2^\i R + h.c. 

~(h* b b R ^tr 2 + h*t R $l) (htiT2$y>* R - ht$2t* R ), 

2\^it 2 Q\ 2 + 2\& 2 Q\ 2 - Q f Q + $£$ 2 ) 



9l 



+ 9 -(H$2 - also [\{Q ] Q) - l(t R t R ) + l(b R l R ) 



[2.3) 



(2.4) 



where g w and g' are the usual SU(2) L and U(l)y gauge couplings. Further, the interaction 
Lagrangian of the Higgs bosons to the top and bottom quarks is given by 

- Acrmions = (h b b R Q T §{ + h.C. ) + ( h t t R Q T lT 2 § 2 + h.C. ) , (2.5) 

where h t and hf, are the top and bottom Yukawa couplings, respectively. 

With the help of the Lagrangians ( |2.3|) -( [2~5"1) , we now proceed with the calculation 
of the one-loop effective potential. More explicitly, in the MS scheme, the one-loop CP- 
violating effective potential is determined by 



C 



:4 

JZ7r q=t ,b 1 i=l,2 



rrit. ( In 



Q 2 



;2.6) 



In (|2.6|) , Cy is the tree- level Lagrangian of the MSSM Higgs potential 



+ A 2 ($ t 2 $ 2 ) 2 + A 3 (^$i)(^ 2 ) + A 4 ($ t 1 $ 2 )($ 2 $i), 



(2.7) 



with 



Ai 



-777,1 



H 2 , 



A 5 



/'I 

'2\ 



■m, 



A, 



m 12 



A, 



(2- 



Further, in (|2.6|) , m 2 (with i — t,b) and m 2 fe (with = h, b\,t 2 , b 2 ) denote the eigenvalues 
of the quark and squark mass matrices JVv M. and A4 2 , respectively, which depend on the 
Higgs background fields. Specifically, M.^ M. reads 



M ] M 



\h t \ 2 \<& 2 \ 2 hth* b <$>liT 2 <$> 2 
h*h h <$>\iT 2 <$>* 2 N 2 |$i| 2 



(2.9) 



with eigenvalues 
1 



t(b) 2 



\h b \ 2 \^\ 2 + \h t \ 2 |$ 2 | 2 +(-) V(N 2 |$i| 2 + M 2 |$ 2 



4|/i t | 2 |/i fc | 2 |$ t 1 $ 



(2.10) 



It is easy to see that, for <fif t2 = 0, ( [2.10| ) simplifies to the known expressions: 
and ml = l^l 2 !^] 2 . 



\ht?W: 



2U0I2 



The (4 x 4) squark mass matrix A4 is more complicated. In the weak basis {Q 



(t L , b L ),U = t* R , D = b* R }, M 2 may be cast in the form: 



(M 2 )^ (M 2 ) 5tht \ 



uu* 



UD* 



(2.11) 



V (M 2 )^ (M 2 ) 3U , {M 2 )~ D ~ D , ) 



with 



M 2 Q 1 2 + |/i fe | 2 $i$i + \h t \ 2 m$ 2 l2 - $2$l 



1 



g 2 w - $ 2 $ 2 



(A* 



2\t 



-h t A t ^ir 2 + h t n*§[ir 2 



'QW* 



(M\ Ut = M 2 + \h t \ 2 $ 2 $ 2 + ^ /2 ($ t 1 $ 1 - $ 2 $ : 



(m 2 : 

(M 2 [ 



1 



'UD* 



M 2 h + - ^ ,2 ($ t i$i - ^$2), 

(At 



2\* 



/i f ^ $fir 2 $ 2 . 



'DC/* 



(2.12) 



Here and in the following, we denote by Mq, M 2 and Ml the soft supersymmetry-breaking 
masses of the third generation of squarks. It is rather difficult to express the four eigenvalues 
m 2 k (q k = ti, bi, t 2 , b 2 ) of M 2 in a simple form. However, as we detail in Appendix A. 3, it is 
not necessary to know the analytic form of m 2 k in order to evaluate the Higgs-boson masses 
and mixing angles |3S]] . Of course, for (pf 2 = 0, the field-dependent squark eigenvalues 
simplify to 



m 



ti(t 2 ) 



/2 



M 2 + M 2 + 2\h t \ 2 \$\ 2 + ^ + 9 n ' 2 i- 1 " 2 ' 



(I0?| 2 -I^l 



+ (-) v/[M 2 - M 2 + x t m\ 2 - |0^| 2 )] 2 + 4N 2 \A$ - /i*0?| 2 



bi(b 2 ) 



M 2 +M 2 + 2\h b \ 2 \^ 2 



gl + g' 2 



+(-) VtM 2 - Ml - x b m 2 - |^| 2 )] 2 + 4|^| 2 \Al$ - /i^| 2 



,(2.13) 



s 



¥ 2 ) 



where x t = \{g 2 w - \g' 2 ) and x b = \{gl 

After having set the stage, we now derive the minimization conditions governing the 
ground state of the MSSM one-loop effective potential and then determine the Higgs-boson 
mass matrices. As usual, we consider the following linear expansion of the Higgs doublets 
$1 and $2 around the ground state: 



$1 



V2 



[vi + 4>i + iai) 



3 « 



V2 



(v 2 + 4>2 + ia, 2 ) 



(2.14) 



where v i and v 2 are the moduli of the vacuum expectation values (VEVs) of the Higgs 
doublets and £ is their relative phase. Following |lQl]lI] , we require the vanishing of the 
total tadpole contributions 
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- 1 



^ai(a 2 ) 



<9£ 
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167T 2 



E E 
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dm 2 

901(2) 



9m 2 

Hi 



m„ 



(2.15) 



9a x 



(2) 



m- 



x I In 



2 

Q 2 



- l 



(2.16) 



and the tadpole derivatives {dm 2 /d(j)i( 2 )), (^ m g fc /^0i(2)) and 



where (m 2 ) 

(dm 2 /dai(2)) are given in ( |A.1|) and ( |A.3|) . Moreover, from ( |A.3|) , we readily see that 
T ai = — tan/3T a2 , with tan/3 = v 2 /v\. This last fact []T0] , |TTjl allows us to perform an 



orthogonal rotation in the space spanned by the 'CP-odd' scalars a\ and a 2 , 

cos (5 — sin f3 
sin p cos f3 





(2.17) 



The Higgs potential then has a flat direction with respect to the G° field, i.e., (dCy / 0G°) = 
0, and the G° field becomes the true would-be Goldstone boson eaten by the longitudinal 
component of the Z boson. 

We observe from ( |2.16|) that a relative phase £ between the two Higgs vacuum ex- 
pectation values is induced radiatively in the MS scheme |I||Il]]. However, we should 
stress that the phase £ is renormalization-scheme dependent. For example, one may adopt 
a renormalization scheme, slightly different from the MS one, in which £ is set to zero 
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order-by-order in perturbation theory JT0[. This can be achieved by requiring the bilinear 
Higgs-mixing mass m\ 2 to be real at the tree level, but to receive an imaginary counter-term 
(CT),Im TTI121 at higher orders, which is determined by the vanishing of the CP-odd tadpole 
parameters T ai and T a2 for £ = 0. As we detail below, the scheme of renormalization of 
Im(m 2 2 e^) does not directly affect the renormalization scheme of other physical kinematic 
parameters of the theory to one loop, such as Higgs-boson masses and tan/3. In fact, it 
has been explicitly demonstrated in QTOJ that physical CP-violating transition amplitudes, 
such as scalar-pseudoscalar transitions, are independent of the renormalization subtraction 
point Q 2 and the choice of phase £. In the following, we adopt the £ = scheme of renor- 
malization, as irrelevant ^-dependent phases in the effective chargino and neutralino mass 
matrices can thereby be completely avoided. 

In the remaining part of this section, we evaluate the one-loop effective Higgs-boson 
mass matrices. Employing the tadpole conditions Tfa = T^ 2 = and T ai = T a2 = allows 
one to substitute the mass parameters [x\ and and lm.m\ 2 into the effective potential 
Q2.6| ). After performing the above substitutions, we can express the charged- Higgs-boson 
mass matrix as follows: 



(M 2 ± h 



5 



Vi 



-2 



V1V2 

ViVj 

dm\ h 



d 2 m 2 q 



Rem? 



.12 



1 9 w viv 2 



i(l - Sij) 



dm 2 

Hk 



dctj 



Vi 



dm 2 



16tt 5 



Qk 



E E 

q=t,b 1 fe=l,2 
2 



171% 

In- qk 



Q 2 



m 2 



- 1 



d 2 m 2 

Qk 



(2.18) 



Since detA^ = in ( |2.18|) , the square of the charged Higgs-boson mass Mjj + may be 
determined by the matrix element {M. 2 ±)i2: 



M 2 



H+ 



VlV 2 



Rem 12 + -g w viv 2 



16vr 2 



q=t,b 



E E 



> / m 5 
Xm <ik[ ln -^t 



d 2 m 2 



fc=l,2 



2 / 1 m q 



d 2 m 2 

Hk 



dm\ k 



v\ \ da,2 



(2.19) 



The self-energy derivatives appearing in (|2.18| ) and (|2.19| ): (d^m^/d^idtp^) and 
(d 2 m 2 /d<pid(p2 ), as well as the tadpole terms (dm 2 k /d(j)j), (dm 2 k /d(f)j) and (dm 2 k /daj), 
are exhibited in Appendix A. 

By analogy, in the weak basis {0i, 02, a i, ^2}, the neutral-Higgs-boson mass matrix 
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takes on the formf] 



Ml 



Ml 



M 2 SP 
Ml 



(M 2 P f Ml 1 ' (2 ' 20) 

where M 2 S , M P and M 2 SP denote the two- by-two matrices of the scalar, pseudoscalar and 
scalar-pseudoscalar squared mass terms of the neutral Higgs bosons. Observe that the 
presence of CP-violating self-energy terms leads to mass eigenstates with no well-defined 
CP quantum numbers. Therefore the CP-odd Higgs-boson mass cannot be identified 
with any of the neutral Higgs-boson masses. The individual matrix elements of Mq are 
given by 
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(2.21] 
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(2.22) 
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dm\ k 



(2.23) 



Again, the analytic expressions for the self-energy and tadpole derivatives with respect to 
the background Higgs fields are given in Appendix A. 

Since G° does not mix with the other neutral fields, the (4 x 4) matrix Mq reduces 
to a (3 x 3) matrix, which we denote by M 2 N . In the weak basis {</>i, </> 2 , a}, the reduced 

^Notice that our convention differs from that given in |Ll| , as the neutral-Higgs-boson mass matrix A4$ 
in that work is expressed in the weak basis {eti, a%, <fii, ^2}- 
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neutral mass-squared matrix M. 2 N may be expressed by 



Ml 



( 



(•M|)n 
(M 2 



S)21 



(M 2 S ) 12 
(M 2 S ) 2 2 



(M 2 SP )i2 

— (M 2 



sin P cos /3 

(•M|)l2 / 



sin {3 
1 



(2.24) 



In writing M. 2 N in ( |2.24|) , we have used the properties of the matrix elements of -M| P : 
(M 2 SP )n = — tan/3 (M 2 SP )i2 and (-M| P )22 = — cot (3 (M 2 SP )2i, and likewise for M 2 P . 

Using the expressions ( |2.19| ) and fl2.24| ), we determine the analytic forms of the RG- 
improved charged and neutral Higgs-boson masses in the next Section. 



3 RG-Improved Higgs-Boson Mass Matrices 

In this Section, we perform a one-loop RG improvement of the squared charged Higgs- 
boson mass Mjj + and of the squared neutral Higgs-boson mass matrix Ai 2 N . The RG 
improvement incorporates all leading two-loop logarithmic corrections to the Higgs-boson 
mass-matrix elements, which were already found in the CP-conserving case to give rise 
to significant contributions to the Higgs-boson masses and couplings. In particular, the 
upper bound on the lightest CP-even Higgs mass was found to be strongly affected by 



the two-loop logarithmic corrections []T~5| , |T6| , p~T| , |T^1 . In carrying out the RG improvement 



we follow the procedure outlined in jT8]|, in which the improvement of the Higgs-boson 
mass-matrix elements was performed by carefully applying the process of decoupling of the 
third-generation squarks. 

Within the framework of the RG approach, the dominant contributions to the Higgs- 
boson mass matrix Ai 2 may be written conceptually as a sum of two terms: 

M\m t ) = M 2 (m t ) + M 2 ' th (m t ). (3.1) 

The first term, A4 2 (m t ), contains the genuine logarithmic contributions which determine 
the whole scale dependence of the one-loop effective potential. These contributions would 
be present, even if the left-right mixing of the stop and sbottom states were absent. The 
second term, A4 2 ' th (mt), describes the threshold effect of the decoupling of the heavier 
stop and sbottom squarks and their respective mixing with the lighter states. At the one- 



loop level, the second term is manifestly scale independent. In (|3 . 1| ) , we are interested in 
evaluating the effective potential at m t , since it has been shown Jl6j that this is the scale 
at which two-loop corrections are minimized. As we explain below, the renormalization of 
the above two contributions must proceed in different ways. 
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Let us denote by Qtb the scale of the heaviest third-generation squark, which we 
assume to be higher than the electroweak scale. In the language of the RG approach, we 
have first to consider that the aforementioned threshold contribution is 'frozen' at the scale 
Qtb, -M 2 ' th (Qtb) = AA4 2 (Qtb), with all the involved kinematic parameters defined at this 
particular scale. Then, we have to rescale the threshold contribution with the anomalous 
dimension factors of the relevant Higgs fields: 

Mlf h (m t ) = /\M%(Q tb ){i;\m t )ij\m t ), (3.2) 

where £,i(m t ) is the anomalous dimension factor of the Hi state to be determined below. 
The one-loop matrix elements AAi 2 j(Q t b) depend on the running quark masses at the scale 
Qtb, which have to be conveniently re-expressed as functions of the corresponding running 
masses at m t . Thus, the anomalous- dimension factors in combination with the one-loop 
relation between the quark masses at scales Q t b and m t yield sizeable two-loop corrections 
to the mass-matrix elements originating from the one-loop threshold effects. 

As was already mentioned, the contribution M. (m t ) of the third-generation squarks 
to the effective potential describes the genuine one- and two-loop leading-logarithmic run- 
ning of the Higgs quartic couplings. In this context, there are two important technical 
details that should be mentioned. First, we notice that, in the MSSM, the tree-level Higgs 
quartic couplings A«, with i = 1,2,3,4, are all proportional to the squared gauge cou- 
plings and g' 2 (cf. (|2.8|)). However, the one- loop (3 functions of \ can generally have 
appreciable values, as they are proportional to the fourth power of the top- and bottom- 
quark Yukawa couplings. As a result, the low-energy values of \ differ significantly from 
their tree-level ones. The RG-improved approach followed here is crucial for implementing 
properly the potentially large logarithmic corrections to the Higgs quartic couplings. 

The second technical remark pertains to the RG evolution of the Higgs quartic cou- 
plings Aj, with i = 5,6,7 (for the notation, see fT^JTll ), which are absent in the Born 



approximation to the MSSM Higgs potential. On field-theoretic grounds, these quartic 
couplings must have vanishing one-loop (3 functions, and cannot be generated by RG run- 
ning. However, these quartic couplings are radiatively induced by threshold effects, and 
have already been taken into account in AA4 2 (Q tb ), given by (|3.2j). 

Following the above discussion, we now proceed with the RG improvement of the 
Higgs-boson mass-matrix elements. To this end, we first need to compute the one-loop 
values of the quartic couplings Aj, with i = 1, 2, 3, 4, where the decoupling of the stop and 
sbottom contributions at their appropriate thresholds is properly taken into account. The 
best way to calculate the latter effects is to consider the logarithmic part of the effective 
potential (|2.6|) in the limit where the squark mixing parameters vanish, i.e., fi = A t = 
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A b = 0. The pertinent one-loop running quartic couplings, denoted by Aj- , may then be 
obtained by 



(i) 



A 



(i) 



A« 



327T 2 

9 



4_ A 2 
4 12/ 



M 2 + m 2 



In 



M 2 + m 2 



Q 2 



32vr 2 



Q 2 

(n 2 - 



y2 \ 2 



- (N - y 



In 



o 2 o' 2 ^ 2 

1^12 _ ^£ _ ^_ 

11 4 12 
M 6 2 + m 2 b 



In 



/M 2 +m 2 



V Q 2 



Q 2 



(3.3) 



9l , 9_ 
4 12 



/2\ 2 



In 



- ( N 2 " y 



# ,2 \ 2 /M 2 + mf 



In 



167T 2 



hi 



^ 2 

M 2 +m 2 



/Mj+m 2 . 
V Q 2 

36 



4 12 



'2 X 2 



111 



/M 2 +m 2 



V Q 2 



M b 2 + m 2 b 

Q 2 



(3.4) 



/2 



12 

12 



fj 



Q 2 

2 



, /max (M? + m, 2 , Ml + m?} 
+ ln( Q l * 



4 



</ 2 ^ 



5 



/2 



9 



12 



12 V 4 



/2 



/2 



4 



M 2 + m 2 



12 V 4 



12 

12 



9 



12 



In 



In 



6 



167T 2 



\h t \ 2 \h b 



111 



+ In 



In 



Q 2 

M 2 Q + m\ 

Q 2 

M 2 Q + mj^ ,j 

Q 2 



M 2 Q + m 2 t 

_ Q 2 

( M 2 +m 2 b 
\ Q 2 _ 
' M/ 2 + m 2 



6 



Q 2 



, (3.5) 



/max (M 2 + m 2 , M 2 + m 2 ) 



9 



Q 2 

2 /M 2 +m 



In 



V Q 2 



. (3.6) 



Moreover, we need to know the one-loop running of the soft supersymmetry-breaking pa- 
rameter Rem 2 2 . Gathering the relevant logarithmic terms present in the effective potential 
rBl), we find 



Re m 2 ^ 



\h t \ 2 Re{^A t ) In 



/max (Mq + ml , M 4 2 + mf 



167T 2 

+ \h b \ 2 Re(fxA b ) In 



/max (Afg + m 2 . , M b 2 + m 2 b ) 



Q 2 



(3.7) 



The analytic form of AAi 2 (Qtb) in Q3.2j) can now be obtained by subtracting the one- 
loop Born-improved mass matrix Ai 2 ^ from its total one-loop contribution Ai 2 ^. Here, 

2 (1) 2 ( 1) 

we have in mind the charged and neutral Higgs-boson mass matrices Ai± and Ai^ 
calculated in Section 2. More explicitly, AAA 2 (Qtb) is given by 



AM 2 (Q tb ) = M 2{1 \Q tb ) - M m Rer4P{Q tb ),\V(Q tb ) 



(3. 
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where Ai 2 ^ represents the tree-level functional form of Ai 2 , expressed in terms of 

2(1) 

and Rem 12 . Furthermore, it is essential to stress again that the kinematic parameters 
involved in (|3.8|) , such as masses and couplings, are evaluated at the scale Q t b- 

Another important ingredient in the RG improvement of the Higgs-boson mass ma- 
trices is the analytic two-loop result for the Higgs quartic couplings Ai,...,A 4 . As has 
been done in (|3.3|) - (|3.6|) , we have to include two-loop leading logarithms, by appropriately 
considering the stop and sbottom thresholds. These two-loop leading logarithmic contribu- 

(2) 

tions to the Higgs quartic couplings, which we denote by A- , can be determined by solving 
iteratively the RG equations |L8| . In this way, we obtain 



A* 



(2) 



6\h h 



,(2) 



X (2) 



(327T 

Q\ht 



2\2 



2\2 



(327T 

Z\h t \ 2 \h b 
(16vr 



h t 2 + 



8gl 



In 2 



In 2 



M 2 Q + m 2 

_Q 2 

M 2 Q + m 2 t 

Q 2 



hr 



+ In 2 



(M% + m 



\^Q 2 
M 2 + m\ 

Q 2 



> (3-9) 
(3.10) 



2^2 



h t \ 2 + \h\ 2 - 8g 2 s 



x 



In 2 



M 2 Q + m 2 

Q 2 



+ In 



2 /max (M t 2 + m 2 t , M b 2 + ml 



Q 2 



(3.11) 
(3.12) 



For later convenience, we define collectively the sum of the tree, one-loop and two-loop 
quartic couplings as follows: 

A, = A, + A, (1) + Xf , (3.13) 

with i = 1,2,3,4. Similarly, the sum of the tree, one-loop and two-loop contributions to 
the soft-bilinear Higgs mixing may be defined as 



Rem 2 2 



Rem 2 2 + Remi9 + Rem 



2(2) 
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(3.14) 



2(2) 

As we see below, knowledge of the two- loop contribution Re m 12 is not required in the 
one-loop RG improvement of the MSSM Higgs potential. 

Given the above definitions of the quartic couplings and the soft Higgs-mixing pa- 
rameter in (|3.13|) and (|3.14j) , we can express the one- and two-loop leading logarithmic 
contributions M. (m t ) to Ai 2 (m t ) by means of the two-loop Born-improved mass matrix: 



M (mt) 



M 2(0) 



Rem 2 2 ( 



m t ),\i(m t ) 



(3.15) 



Note that M. (mt) also includes the tree-level terms. As has also been stated explicitly in 
( |3.15| ), M. (m t ) is expressed in terms of mass and coupling parameters evaluated at the 
top-quark-mass scale. 
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The last ingredient for completing the programme of RG improvement of the Higgs- 
boson mass matrices is knowledge of the analytic expressions for the anomalous dimension 
factors that occur in ( |3.2| ). These analytic expressions are given below for the charged and 
neutral Higgs-boson cases separately 

Adopting the framework outlined above []TS| , it is not difficult to compute the RG- 
improved charged Higgs-boson mass M^ + (m t ) at the top-mass scale through the relation: 

M 2 H+ (m t ) = M 2 H+ (m t ) + [£+(m t ) £ 2 ~(m t )] (AM 2 H+ f{Q tb ) + (M^+) th (rrit) , (3.16) 

where Q 2 tb = max (^Mq+tti 2 , M 2 +m 2 , M 2 +m^j, and ^ + (m t ) and i 2 {m t ) are the anomalous 
dimension factors of the charged Higgs fields (ftf and <p 2 , respectively: 



?1+(roi ). 1 + |^ ln |, &(mi) . 1 + |^ hl |. (3 ,T) 



Further, M H+ (m t ) is the squared two-loop Born-improved charged Higgs-boson mass given 
by 

and (AM 2 I+ ) tb is the one-loop scale-invariant part that contains the stop and sbottom 
contributions: 

T 2 \ib 



(AM 2 H+ nQ tb ) = - 1 A {M^ 1] t 2 {Q th ) + Rem 2 £\Q tb) 

sm p(m t ) cos p{m t ) 1 

1 \0), 



+ -\r(Q tb )v l (Q tb )v 2 (Q tb )\, (3.19) 

where the scale at which the kinematic parameters are to be evaluated has been indicated 
explicitly. In ( |3.19| ), the term between brackets is the threshold contribution to the off- 



diagonal matrix element of the charged-Higgs-boson mass matrix, where (-M±^)^| denotes 
the one-loop contribution of the third-generation squarks to (M±)i 2 . Finally, {M 2 I j+ > ) th 
describes the one- loop quark contribution to (see ( |3.16| )). 

We remark that the charged-Higgs-boson mass matrix receives the common anoma- 
lous dimension factor "(m t )^ "( m t) , even though different matrix elements of M. 2 ± are 
involved. This is because hA\ must possess a vanishing determinant at any RG scale 
Q 2 , as one of its mass eigenstates must correspond to the massless would-be Goldstone 
boson G + that forms the longitudinal component of the W + boson. As a consequence, 
the following relations among the matrix elements of the RG-frozen part AM. 2 ± of the 
charged-Higgs-boson mass matrix are obtained: 

(AM 2 ± ) n (Q tb ) = - tan /3(Q tb ) (AM 2 ± ) 12 (Qtb) , 

(AM 2 ± ) 22 (Q tb ) = -cotp(Q tb ) (AMl) 21 {Q tb ) . (3.20) 
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After including the RG running due to the Higgs-boson anomalous dimensions, we find 

: — tan/3(m t ) 



(AM 2 ± )n(Qtb) 



(AMl) 12 (Q tb ) 



[£~K)P 

where we have used the RG relation: 

tan^(Q t6 ) 



cot f3{m t ) 



(AMj) 12 (Qtb) 



tan f3{m t ) . 



(3.21] 



(3.22) 



As a consequence of this last relation, it is evident that the RG running of the different 
matrix elements of AAi± may be expressed in terms of the running of (AM 2 ^)^ and the 
value of tan/? at the scale m t . 

Correspondingly, the RG- improved neutral-Higgs-boson mass matrix A4 2 N may be 
computed by 



(M 2 N Um t ) + {AM 2 N %{Q t ) + [4(m t ) 

+ (-M^ {1) )8K) , 



(AM 2 N )l(Q b ) 
(3.23) 



with i, j = 1, 2, 3, Q\ = max + m 2 , M 2 + mfj and Q 2 b = max + m\, Ml + m\ 
Notice that, unlike the charged-Higgs-boson case, one has to introduce here two decoupling 
scales Q t and Qb, as the stop/top and sbottom/bottom loop effects occur separately in the 
threshold contributions. The parameters ^ (with q = t,b) are the anomalous- dimension 
factors related to the neutral Higgs-boson fields 



&( m t) £j(™t) , for 
€i( m t) £>Oi) > fo r 



i,j = 1,2 

i = 1, 2, 3 and j = 3 



(3.24) 
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1 + 



1 + 



In 



Ql 



327T 2 

3\h b \ 2 ,_ Ql 



2 ' 



32tt 2 



In 



2 ' 



3\h t \ 2 Q 2 



In 



m 



1 + 



32tt 2 



3N 2 Q 2 b 



2 ' 



327T 2 



In 



mi 



(3.25) 



We should observe that, for the very same reasons as in the charged-Higgs-boson case, 
the vanishing of the determinants of M. 2 P and M. 2 SP at any Q 2 scale leads to the com- 
mon anomalous dimension factor ^(mt) = ^i{m t ) = ^(^t) in the calculation of 



Ai 2 N (m t ) in fl3.23|) . Since this last fact involves the matrix elements (A4 2 P )i 2 , (A / i| P )i2 and 
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(M 2 SP )2i 1 the corresponding matrix elements of (AvM|r)f 3 in ( f3.23| ) are given by 

cos p(m t ) 

(AM%%(Q q ) = --^^MizTT, (3-26) 



(M 2 SP )UQ q ) / A \9 tr, n (-MIp)Ii(^) 



sin f3(m t ) 



sin /3(mt) cos f3{m-t) ' 

where the RG-scale dependence of the involved quantities has been displayed explicitly. 

In (|3.23|) , the (3 x 3) matrices A4 N and (AAI^)' (with q = t,b) describe respectively 
the two-loop Born-improved effects and the one-loop threshold contributions associated 
with the decoupling of the heavy squark states: 



M N (m t ) 
(AM 2 N y(Q q ) 



M 



2(0) 
JV 



Remf 2 (m t ), \i{m t ), A 2 (m t ), A 34 (m 4 ) 



(3.27) 



(M 2 N W Y(Q q ) ~ MT 



where A34 = A 3 ±j + A 4 1; (likewise A34 = A3 + A4) and A4^ v> is the tree-level functional form 
of M 2 ■ 



Reml^(Q q ), X?'%Q q ), A^Q,), A&^Q,) ] ,(3-28) 



2(°) ; r 



N' 



Remf 2 
sin (3 cos (3 



\ 



sin 2 (3 — sin (3 cos (3 \ 
sin /3 cos (3 cos 2 /? 
1/ 



— v 



(3.29) 

/ 2AiCos 2 /9 A 34 sin /3 cos (3 \ 

A34 sin (3 cos /? 2A 2 sin 2 (3 , 
v / 



with A34 = A3 + A 4 . As in the charged-Higgs-boson case, we write {M.^ ) q to denote 
the one- loop part of Ai% containing the contributions of the third-generation squarks, and 
(.M 2 ^)^ to denote its fermionic one-loop counterpart. 

The resulting RG-improved Higgs-boson mass matrix M. 2 N (m t ) is a symmetric, 
positive-definite (3 x 3) matrix, and can therefore be diagonalized by an orthogonal trans- 
formation as follows: 



T M 2 N {m t )0 = diag M 2 (mt), M 2 (m t ), M 2 {m t ) 



(3.30) 



where we have defined the Higgs fields such that their RG-improved masses satisfy the 
inequality: 

M Hl {m t ) < M H2 {m t ) < M Ha (m t ) . (3.31) 

Notice that our convention in ( |3.30|) differs from that chosen in |Tl] , as we assign the Higgs 
fields in the reversed order. Analytic expressions for M^(m t ) and O are presented in 
Appendix B. 
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Before closing this Section, two important remarks are in order. First, we observe 
that the free kinematic parameters of the MSSM Higgs sector are 



M H+ (m t ) 
M 2 Q (Qt 



fib 



tan(3(m t ), n(Qtb 
M 2 (Q tb ), M 2 (Q 



At(Qtb) , A b (Qtb) 



tb) 



(3.32) 



In fact, the soft Higgs-mixing parameter Rem 2 2 (m t ) may be substituted by the squared 
RG-improved mass Mjj+(m t ) of the charged Higgs boson (cf. ( BjT§ ) and in the 



neutral Higgs-boson mass matrix Ai%(m t ) in ( 3.23|) 



Secondly, we reiterate the fact that Im m\ 2 can be renormalized independently, with- 



out affecting the renormalization of the physical parameters of the theory fllQf . As was 
stressed in Section 3, the £ = scheme of renormalization gives rise to a considerable sim- 
plification, since we can get rid of the radiatively-induced phase £ between the two Higgs 
vacuum expectation values in the analytic expressions of the Higgs-boson masses and mix- 
ing angles. For example, within the above £ = scheme, the mass renormalization of H + 
may be entirely reabsorbed by a corresponding renormalization of Rem 2 2 and A4. In other 
words, it can be shown that Mjj+ is Q 2 -independent, after including the RG running of 
Rc 

7TI12 an d A4, denoted as 7Rem 2 an d Pxt- For simplicity, we assume that only the third 
generation of squarks contributes to 7R em 2 2 , since fermions do not contribute to the RG 
running of Imm^. The analytic forms of 7R e mf 2 an d /3a 4 axe given by 

3 



7Re 



'12 



16vr 2 



/5a 4 = 
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\h t \ 2 Re({iA t 
3 



I hb I Re (fiAt 



rflnQ 2 



16tt 2 



2\h t \ 2 \h b \ 



9 j(\h 



, 2 -1- 
t + 



h h \ 2 



+ 



9t 



(3.33) 
(3.34) 



Obviously, the RG running of Remf 2 due to t and b is only relevant for non-zero values 



of fx At and iiA b . Employing ( p. 19 ) and examining only the In Q 2 -dependent part, one can 
verify that 

'8 2 TrM i \ 1 /dTr.M 4 ' 



dM 2 H+ 



0, (3.35) 



d\nQ 2 ! ^ em ^ 2' A4 1 ' 32vr 2 V\ d^t dfa I v x \ da 2 

as it should be. As can also be seen from (|3.35|) , an important role in this proof is played 
by the necessary CP-odd tadpole term (<9Tr M. A jda-i) [0. 



4 Effective Top and Bottom Yukawa Couplings 

In addition to the RG improvement of the Higgs-boson mass matrices discussed in the 
previous Section, we consider here a further improvement related to the non-logarithmic 
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threshold corrections to the top- and bottom-quark Yukawa couplings. Specifically, apart 
from the usual RG running, the effective top- and bottom-quark Yukawa couplings obtain 
additional non-logarithmic threshold contributions, which are induced by the decoupling 
of the heavy SUSY states at a high scale, e.g., Qtb- For the bottom-quark Yukawa case, 
the one-loop RG relation between the bottom mass and the bottom-quark Yukawa cou- 
pling at the scale Q t b receives quantum corrections that also include terms proportional to 
tan (3 . Since these last terms can be significant for large values of tan (3^ we must 

resum them within the RG approach, so that the actual size of the radiative corrections to 
the Higgs-boson masses and couplings can properly be extracted [l9|,|42j . For the top-quark 
case, instead, although one-loop suppressed, the respective corrections can still give rise to 



an enhancement of up to 4 GeV in the prediction for the lightest Higgs-boson mass |43|J44 
and therefore should be included in the computation. 

There may also be important CP-violating one-loop corrections to the bottom- and 
top-quark Yukawa couplings, in addition to the CP-violating effects induced by the radiative 
mixing of the Higgs states, which were considered in Sections 2 and 3 in detail. In the 
leptonic sector, these CP-violating vertex corrections are generally small However, 
the CP-violating radiative corrections to the couplings of the Higgs bosons to b quarks are 
significant [fCll , because of the large Yukawa and colour-enhanced QCD interactions [39 



In particular, the radiative effects of the Higgs-boson couplings to the bottom quarks can 
be further enhanced, if the respective Higgs-mass eigenstate couples predominantly to the 
Higgs doublet $ 2 [IIIED], as the tree-level 6-quark Yukawa coupling is suppressed in this 
case. For a general discussion of the form and the origin of these finite Yukawa corrections 
to the third-generation quark masses, the reader is referred to the original literature |40|j4T 



In the following, we give a brief discussion of the non-logarithmic corrections to the top and 
bottom Yukawa couplings, and pay special attention to the CP-violating vertex effects. 

We start our discussion by considering the effective Lagrangian of the 6-quark Yukawa 
coupling [|H||TT||: 



i>?b R b L 



Ah b b R b L 



h.c. , 



(4.1) 



with 



5hb 
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h h 



2c^s * a j, 2 2 I I 

m~A b l(mr , ml , m 5 



2a, 



m*^*I{ml x1 m 



2 I |2\ 
b 2 > \ m ~9\ ) 



\2jl 2 2 I I 

H\ I{m k ,m h , \n\ 



167T 2 
17) I 2 

\' b t\ a* *rf 2 2 I |2\ 

I{m~ ti ,m i2 ,\n\ ) 



(4.2) 
(4.3) 



§For the r-lepton Yukawa coupling, the corresponding enhanced tan/3 terms are much smaller, because 
they are proportional to the weak gauge couplings |t(],[l^] . 
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where ot s = g 2 /(4:7i) is the SU(3) C coupling strength, and I (a, b, c) is the one-loop function 

abln(a/b) + bc\n(b/c) + ac\ia(c/a) . . 

(a — b)(b — c)(a — c) 

The 6-quark Yukawa coupling h b (Q tb ) is then related to the running 6-quark mass m b (Q tb ) 
by 

b ' ^/2M W cos f3 [I + 5h b /h b + (Ah b /h b ) tan p] ' 1 ' ' 

where Sh b /h b and Ah b /h b are given in (|4.2j ) and (|4.3|) , respectively. The running 6-quark 
mass m b (Q tb ) is obtained by means of the RG running of the 6-quark mass from the scale m b . 
In ( |4.5| ), we have redefined the right-handed 6-quark superfield, so that the physical 6-quark 
mass is positive. Under such a field redefinition, only the Yukawa coupling h b becomes com- 
plex, while the phases of 5h b /h b and Ah b /h b as well as those of the supersymmetry-breaking 
parameters do not change. Moreover, since only the moduli of the Yukawa couplings h b 
and ht enter the field-dependent quark and squark masses in (|2.10|) and (|2.13|) , the neutral 
Higgs-boson mass matrices remain unaffected by the above field redefinition. Also, we have 
checked that the very same property of invariance under rephasings of h t and h b persists for 
the charged Higgs-boson mass matrix as well. At this point, it is interesting to observe that 
the sum 5h b + Ah b tan /3 in ( |4.5| ) receives two sorts of quantum corrections, one originating 
from QCD effects and another from a chargino-mediated graph. The QCD correction is pro- 
portional to the hermitean conjugate of the sbottom-mixing parameter X b = A b — fi* tan j3, 
whilst the char gino- induced diagram |3D| depends linearly on the stop-mixing parameter 
X t = A t - n* cot (3. 

The effective Lagrangian describing the t-quark Yukawa coupling is given by 

- Z> ft = Ah t 4>lt R t L + {h t + 8ht)4>lt R t L + h.c. (4.6) 

The corresponding relation for ht as a function of m t may easily be determined analogously 
by the effective Lagrangian (O), and reads 



h _ 9wm t , 4? , 

* ' V2M W sin (3 [1 + 5h t /h t + (Ah t /h t ) cot (3} ' 1 ' ' 



with 
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As in the case of the 6-quark Yukawa coupling, we have to make a judicious phase rotation 
of the right-handed t-quark superfield, such that the physical top-quark mass becomes 
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positive. Again, one can show that such a field redefinition does not change the analytic 
results of the RG analysis. 

At this stage, it is important to remark that, within the RG-resummation approach 
described in Section 3, the non-logarithmic corrections must be treated as threshold effects 
and hence they should only contribute to the RG-frozen part of the Higgs-boson mass 
matrices, generically denoted as AM 2 (Qtb)- Therefore, the decoupling procedure for the 
heavy squark states requires that the effective b- and t-quark Yukawa couplings given by 
( |15D and (|4.7|) are evaluated at the scale Q t b- As we discuss in Section 5, these additional 
Yukawa corrections can lead to observable effects in Higgs-boson searches. 

It is now straightforward to obtain the interaction Lagrangians of the Higgs-boson 
mass eigenstates Hi to the up- and down-type quarks, collectively denoted as u and d. 
Taking into account both CP-violating self-energy and vertex effects, we find 
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(4.10) 



(4.11) 
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(4.14) 



where the Higgs scalar and pseudoscalar couplings are normalized with respect to their SM 
values. 

Finally, it is interesting to investigate the behaviour of self-energy- and vertex-type 
CP violation in the decoupling limit of a heavy charged Higgs boson in the MSSM. Thus, 
for values of the charged Higgs mass M H + 3> Mz, one has O31 — > 0, while On — > cos (3 and 
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O21 - ► sin (3. In this limit, the scalar components of the H\dd and H\uu couplings acquire 
the known SM form, given by g w m d / (2M W ) and g w m u /(2M w ), respectively, where 

m d = -j= (h d + 5h d + A/i d tan/5) v 1 , 

m u = (h u + 5h u + Ah u cot (3 ) v 2 (4-15) 

have already been defined to be positive in ( |4.5|) and ( |4.7|) . For similar reasons, the pseu- 
doscalar parts of the H\dd and H\uu couplings vanish, as they are proportional to O31 and 
lmm du = 0. On the other hand, in the same large M#+ limit, the scalar and pseudoscalar 
couplings of both the two heaviest Higgs bosons H 2 and H 3 to the up and down fermions 
do not vanish. We can therefore conclude that CP-violating self-energy and vertex effects 
do not decouple in the heavy neutral Higgs sector. In the next section, we demonstrate ex- 
plicitly the aforementioned (non-) decoupling features of CP violation by analyzing specific 
phenomenological examples. 



5 Phenomenological Discussion 

In this Section we discuss the phenomenological implications of radiative Higgs-sector CP 
violation in the MSSM for Higgs-boson searches at high-energy colliders. We focus our 
attention on the physics potential for discovering Higgs bosons with mixed CP parities 
at LEP 2 and the upgraded Tevatron collider, and also comment on the enhanced search 
capabilities offered by the LHC. 

At the LEP 2 and Tevatron colliders, neutral Higgs bosons are predominantly pro- 
duced via the Higgs-strahlung processes in e + e~ and qq collisions, such as e + e~ — ► Z* — > 
ZHi g§, qq -> Z* -> ZHi and qq -> H /± * -> W ± H i g7J, with i = 1, 2, 3. If the next-to- 
lightest Higgs boson is not too heavy, Higgs bosons can also be produced copiously in pairs 
through the reactions: e + e~ — > Z* — » HiHj and qq —>■ Z* — > HiHj. In addition to the 
Higgs-boson masses, the Higgs-boson couplings to the gauge fields play an essential role 
in our forthcoming discussion. The effective Lagrangians governing the interactions of the 
Higgs bosons with the W ± and Z bosons are given by 



C HVV = g w M w j2 9H i vv(H i W+W-'> 1 + - \ H^Z") , (5.1) 



i=1 2 cos 2 9 U 

3 

Chhtw± = ^J2 9H iH -w + (Hii\H-)W + ^ + h.c, (5.2) 

i=l 
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where cos 6^ = M w /M z , = — d^, and 

QHiW = cos (3 O u + smf30 2i , (5.4) 

9HiHjZ = 3i (cos P 2j - sin flOy) - 3j (cos/30 2 i - sin^OijJ , (5.5) 

9HiH-w+ = cosp0 2i - sin (3 On + i0 3i . (5.6) 



For completeness, we have included in ( ]5.2[ ) the interactions of the charged Higgs bosons 
H ± with the neutral Higgs and W T bosons. Note that the couplings HiZZ and HiW + W~ 
are related to the HiHjZ couplings through 

9H k VV = ZijkdHiHjZ- (5.7) 

Moreover, unitarity provides the constraint 

E 9%vv = 1 • (5-8) 
i=i 

Evidently, if two Higgs-boson couplings to gauge bosons are known, this is sufficient to 
determine the complete set of the couplings QHiW and gn.H.z 



In the above calculation of the effective Higgs-gauge-boson couplings, we have 
assumed that the dominant contributions arise from Higgs-mixing effects. As op- 
posed to the 6-quark Yukawa case discussed in Section 4, proper vertex correc- 
tions to the HiZZ and HiHjZ couplings do not contain strong-coupling- or tan/3- 
enhanced diagrams. Therefore, naive dimensional analysis suggests that these cor- 
rections are suppressed relative to their tree-level values by loop factors of the kind: 
(3a m /47r) (1/iAtl/m?) , (3a w / Air) {\iiA t \ 2 v 2 / m®^) < 1CT 2 , where t\ is the heaviest stop 
squark (see Q2.131 )). In the following, we neglect proper vertex corrections to the HiZZ and 
HiHjZ couplings. 

For our phenomenological discussion of CP violation, we consider the following two 
representative values for tan/?: (i) tan/3 = 4 and (ii) tan/3 = 20. For definiteness, unless 
stated otherwise, the soft supersymmetry-breaking and /i parameters are set to the values 

Msusy = M Q = M t = M b = 0.5 TeV, /i = 2TeV, \A t \ = \A b \ = ITeV, 
\rrig\ = \m^\ = 0.3 TeV, \ m g\ — 1 TeV , (5.9) 



As can be seen in ( |5.9| ), we have chosen relatively large values for the stop and sbottom 
mixing parameters A t , A b and /i, as well as a common left- and right-handed squark mass 
^susy, which leads to enhanced CP- violating effects of the CP-odd phases arg (pA t) b) on 
the Higgs sector. 



24 



As was mentioned in the Introduction, large CP-odd phases may lead to rather large 
EDM contributions, thereby violating the known upper bounds on the electron and neutron 
EDMs d e and d n : d e /e < 0.5 10~ 26 cm gjj and d n /e < 0.6 10~ 25 cm [[?| at the 2-a level. 
One phenomenologically interesting possibility for avoiding the possible CP crisis is to make 
the first two generations of squarks rather heavy with masses much above the TeV scale 
27,281, keeping the third generation relatively light with masses of order 0.5 TeV. In such 



a scenario, CP violation may only reside in the third generation. For our illustrations, 
we shall take the [i parameter to be real and assume that the only CP-odd phases in the 
theory are arg (At) = arg (Ab) and arg (frig). Again, as one way to avoid the one-loop EDM 
constraints |29j, we have taken a gluino mass of 1 TeV in ( |5.9| ). However, in such a scheme, 
one has to worry about the fact that Higgs-boson two-loop contributions to the EDMs 
|3"4| might still become sizeable. For the low-tan (3 scenario in ( |5.9| ), the two- loop EDM 
contributions are of the order of the experimental EDM upper bounds mentioned above. 
Since these two-loop EDM effects depend almost linearly on tan/3, one might then need 
to arrange for cancellations among the different one- and two-loop EDM terms |29[ at the 
level of 10% for the scenario with tan (3 = 20. We believe that this can be achieved without 
excessive fine-tuning of the CP-violating parameters of the theory. 

As was already noticed in [ 14, 18, 11 1, the radiative corrections to the lightest Higgs 
boson Hi depend crucially on the stop mixing parameter \X t \ = \A t — p? cot (3\. Specifically, 
the radiatively-corrected ifi-boson mass increases as \X t \ increases, reaching a maximum 
when |X 4 |/Msusy ~ 2.45. Then, as \X t \ further increases, the radiative corrections to 
iJi-boson mass decrease and may even become negative, driving the latter to very small, 
experimentally excluded values. A distinctive feature of the CP-violating SUSY scenario 
compared to the CP-conserving one is that \X t \ can be increased by varying only the phase 
arg (A t ) from zero to higher values, but holding fixed \ A t \ and For similar reasons, high 
values of \X t \ induced by large values of arg (A t ) can make the mass of the lightest stop 
squark t% very low, so as to violate present experimental constraints, i.e., mf 2 Z 100 GeV. 
Furthermore, light stop quarks, with rrif < 300 GeV and large \X t \ values, can give rise 
to observably large contributions to the electroweak precision parameter Ap ||49|| . For the 
scenarios under discussion, mi 2 is always larger than about 300 GeV, for all the parameter 
space for which the i^i-boson mass acquires acceptable values. Therefore, apart from the 
bounds derived by EDM constraints, the requirement that the lightest Higgs-boson mass 
is positive can be used naively to set an upper bound on the phase of A t . 

In Fig. [I] we give numerical predictions for the two lightest Higgs-boson masses 
and M H2 , and for the three relevant HiZZ couplings squared as a function of arg(^4 t ), 
for two different values of &rg(rrig): arg(mg) = (solid lines) and arg(mg) = n/2 (dashed 
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lines). We first discuss the scenario with tan/3 = 4, for which the values of the remaining 
soft supersymmetry-breaking parameters and \x are given in (|5.9| ). Since our interest is to 
analyze dominant CP-violating effects for a light Higgs sector, we present predictions for a 
relatively small charged- Higgs-boson mass, = 150 GeV. In the CP-conserving limit of 
the theory (arg (A t ) = 0), the mass M Hl of the lightest neutral Higgs boson is close to 85 
GeV, while the square of the H X ZZ coupling, g 2 HlZZ) is approximately equal to 0.8. These 



values of masses and couplings are now excluded by Higgs-boson searches at LEP 2 f22 
However, this situation changes crucially once CP- violating phases become relevant. As the 
phase of A t increases, two important effects take place. First, as was mentioned above, the 
stop mixing parameter \X t \ becomes larger, giving rise to larger ifi-boson masses. Second, 
the mass-matrix terms describing the scalar-pseudoscalar mixing are enhanced, thereby 
effectively leading to large modifications in the couplings of the Higgs bosons to gauge 
bosons. This second effect can be attributed entirely to CP violation. In fact, as can be 
seen from Fig. |l](b), for arg (A t ) « 80 degrees, g 2 HlZ z gets YeT J suppressed, implying that 
LEP 2 cannot detect the Higgs boson Hi via e + e~ — ► Z* — > ZHi. On the other hand, 
for the same range of values of org (At), i.e., axg(A t ) = 80°-95°, the H2 and H$ bosons 
have significant couplings to the Z bosons. Although the H% boson is too heavy to be 
detected at LEP 2 in this case, the H 2 boson has a mass of 105-110 GeV and g\ 2 zz ~ 
0.6, which may be probed at LEP 2 in this year's run. For larger values of the phase of 
A t , 95° < arg (A t ) < 110°, the discovery of a Higgs boson at LEP 2 is more challenging. 
The lightest Higgs boson Hi acquires a mass below 90 GeV, but the HiZZ coupling is 
too small to allow experimental detection through the reaction e + e~ — > Z* — > ZH\. In 
addition, the H 2 boson becomes too heavy to allow for discovery via e + e~ — > Z* — > ZH 2 , 

with 9h 2 zz ~ °- 7 - 



As was discussed in |TTJ, the Hi and H 2 bosons may also be searched for in the 
channel e + e~ — > Z* — > HiH 2 . Since the squared coupling g 2 HlH2Z = 9h 3 zz ~ 0-2 almost 
independently of arg (At), a careful experimental analysis of the parameter region of interest 
will be necessary to determine whether the two lightest Higgs bosons can be observed for 
such large mass differences (M Hl — M H2 > 40 GeV) and such small gH 1 H 2 z couplings. 

In Fig. [I] we also present predictions for a gluino phase of 90 degrees (dashed lines). 
Since the vertex corrections are generally small for low or moderate values of tan/?, they 
are expected to induce only small corrections to the Higgs-boson masses and mixings. This 
last fact is reflected in Fig. [I], even though the coupling g 2 HlZ z (.9h 2 zz) S e ^ s slightly smaller 
(larger) for larger values of the phase of A t . 

Fig. ^| shows the changes in the predictions for the same choice of parameters as in 
Fig. [I], but with M SU sy, A* an d \A t \ = \A b \ rescaled by a factor of 2. This rescaling leads 



26 



to a slight increase (decrease) of the Higgs-boson mass Mr 1 (Mh 2 ), while g 2 HlZ z exhibits a 
slightly different quantitative dependence on the phase of A t , especially for the region in 
which g'jj 1 zz is very small. Thus, although Mh 1 becomes small for arg (A t ) > 115 degrees, 
the H\ZZ coupling gets sizeable again, well within the capabilities of LEP 2 to test. 

We now investigate more quantitatively the predictions of a large-tan (3 scenario 
for the Higgs-boson masses and couplings. We adopt the scenario given in (|5.9|), with 
tan/3 = 20 and M H + = 150 GeV. In this large-tan/? scenario, one has |//|cot/3 <C \A t \, 
and the effective stop mixing parameter \X t \ ~ \A t \ is almost independent of a.rg(A t ). 
Therefore, as is seen in Fig. ||](a), the Higgs-boson masses Mh 1 and Mh 2 do not exhibit 
any significant variation as a function of axg(A t ). In contrast to the Higgs-boson masses, 
Fig. |^(b) shows that there is a non-trivial dependence of the squared couplings g 2 H . zz on 
arg(v4f). Furthermore, the next-to-lightest Higgs boson H 2 is heavy enough to render its 
search through the e + e~ — > Z* — > ZH2 reaction kinematically inaccessible at LEP 2. For 
similar reasons, we find that, for all values of axg(A t ), the Hi and H 2 bosons are rather 
too heavy to be produced via the e + e~ — > Z* — > HiH 2 channel at LEP 2. As a re- 
sult, Higgs-boson searches at LEP 2 tend to be more efficient for small values of the A t 
phases, for which the lightest neutral Higgs-boson mass is close to 100 GeV and g 2 HlZZ * s 
non- negligible (g 2 H lZZ ~ 0-3)- In addition, for large A t phases, arg (A t ) ^ 80°, the HiVV 
coupling (with V — Z, W) is rather suppressed, so that the H x Higgs boson, although it 
becomes lighter with a mass in the range 90-95 GeV, will be elusive at LEP 2, and may 
also escape detection via the corresponding channel at the upgraded Tevatron. However, 
the next-to-lightest Higgs boson H 2 has couplings of order unity to the Z and W bosons. 
Present simulations show that a neutral Higgs boson, such as H 2 , with Mh 2 ~ 180 GeV and 
a SM-like coupling strength to vector gauge bosons can be tested at the Tevatron collider 
with a total integrated luminosity of 10 fb -1 . However, discovery of such a Higgs boson 
at the 5-cr level would demand a total integrated luminosity of 30 fb -1 , and would have 
a reach up to Mh 2 ~ 120 GeV [p0| . Finally, even though the HiH 3 Z coupling is close to 
unity (g^HzZ = 9h 2 zz) f° r ar S (At) > 100°, the Higgs-pair production of Hi and H 3 is not 
kinematically allowed at LEP 2, since M^ 3 w M H +. Further studies will be necessary to 
investigate the potential of this production mechanism at the Tevatron. 

It is interesting to present predictions for the neutral Higgs-boson masses and their 
couplings to gauge bosons for lower values of the charged Higgs-boson mass M H + in the 
above large-tan (3 scenario. In Fig. f|, we plot numerical estimates for the same kinematic 
parameters as in Fig. |3], but with Mh+ = 135 GeV. In this case, the ifi-boson mass varies 
approximately between 80 and 65 GeV, and the HiZZ coupling rapidly decreases as the 
phase of A t increases. The two heaviest neutral Higgs bosons H 2 and H 3 have masses in the 
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range between 120 and 130 GeV. Hence, these two Higgs bosons cannot be produced via 
e + e" — > Z* — > ZH 2 or e + e~ — > Z* — > ZH 3 at LEP 2. However, the and H 3 bosons may 
still be accessed via e + e~~ — > Z* — > HiH 2 or HiH 3 . Interestingly, the squared couplings 
9H1H2Z = 9h 3 zz and ^!if 3 z = fi^zz exhibit a cross-over as a function of arg(A t ). The 
crossing point of the two squared couplings is when arg (A t ) ~ 90°. For arg (A t ) = 180°, 
one of the squared couplings goes to and the other to 1, depending on the phase of the 
gluino mass. In this case, the two heaviest neutral Higgs bosons become almost degenerate 
in mass. For the whole range of values of arg (At), either the H 2 or ^3 Higgs boson can be 
tested at the upgraded Tevatron collider provided a total integrated luminosity of 10 ftr 1 
per detector is available p| 

Figs. |^(a) and (b) show the degree of mass splitting between the two heaviest neutral 
Higgs bosons H 2 and if 3, for the same choice of parameters as in Figs. [I] and but for 
M H ± = 200, 300, 400, 500 GeV. As was already observed in (0,0, even though the H 2 
and H 3 bosons are almost degenerate in the CP-conserving limit of the theory, they can 
have a degree of splitting up to 30% for a maximal CP- violating phase arg (A t ) ~ 90°. The 
comparison of the Fig. |5|(a) with (b) reveals that this last result is almost independent of 
the common scale factor of Msusy, /•* and \A t \. Also, the degree of mass splitting is not 
much affected by the value of the gluino phase arg(mg): the results for arg (rrig) = 90° 
are slightly higher than those of arg (nig) = 0. In this vein, it is interesting to mention 
that large CP-violating scalar-pseudoscalar mixings can lead to observable phenomena of 
resonant CP violation at high-energy colliders [^l| ,^2]]. 

In Figs. H and |7|, we examine the behaviours of the scalar and pseudoscalar parts 
of the Hibb coupling as functions of the CP-odd phase axg(A t ), for two different charged 



Higgs-boson masses, Mh+ = 150 and 300 GeV. As was done in |Q|, we find that the 



best way of analyzing such a behaviour is in terms of the CP-even and CP-odd quanti- 
ses: [(g s Hlbb ) 2 + (9H lbb ) 2 } and 2g s Hibb g^ ibb /[(g s Hibb ) 2 + (g% lbb ) 2 ]. For example, Higgs-boson 
branching ratios are proportional to the first quantity, while the second one will only occur 
in CP-violating observables. In other words, ^gfi^bQH^b/^dH^b) 2 + (dn^b) 2 } gives a mea- 
sure of the CP-violating component in the Hibb coupling. If we compare the predictions 
for M H + = 150 GeV with those for M H + = 300 GeV in Figs. | and 0, we find that the 
CP-violating component of the Hibb coupling reduces in magnitude, for large values of the 
charged Higgs-boson mass. Such a decoupling behaviour of the CP-violating Hibb compo- 
nent is in agreement with our observation, which we already made at the end of Section 4. 
From Figs. |6|(b) and 0(b), we see that the impact of the gluino phase on the CP- violating 
component of the Hibb coupling is more important for the large-tan f3 scenario. This may be 
attributed to the fact that the radiatively-induced term Ahb tan ft, which crucially depends 
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on the gluino phase and tan /3, has a dominant contribution to the H\bb coupling. 

There can be a cancellation or a strong suppression of the coupling of the lightest 
Higgs boson Hi to the bottom quarks, depending on the magnitude of the CP-violating 
phases and of the products A t fi, A^n and rrigfi. This cancellation usually takes place for 
moderate values of the charged Higgs mass and large values of tan (5. Such an effect is also 
present in the CP-conserving case, for specific signs and magnitudes of the above products 
involving the trilinear terms A t ^ and the gluino mass, and has been discussed in detail in 
||53|| . Figure |8](b) illustrates such a cancellation for the CP- violating SUSY model under 
discussion. For example, we observe that for the set of SUSY parameters considered in Fig. 
||, the H\bb coupling can be strongly suppressed for arg (A t ) = arg (Aj,) ~ 15°. Moreover, 
we have checked that for this same set of parameters the H\ZZ coupling is almost SM- 
like. In addition, Fig. ||(a) shows that the mass of the lightest Higgs boson is of order 105 
GeV, practically independent of the CP-violating phase. This is therefore an extremely 
interesting example, since the lightest Higgs mass is in the mass range that may be within 
the reach of LEP 2, and its production cross section will be SM-like. However, the main 
decay channel, Hi — > bb, can be strongly suppressed if the CP-violating phases arg(A t ) 
and arg(mp) lie in a specific range. Therefore, in such a scenario, the detection of the Hi 
boson may in principle be impossible, even in the final run of LEP 2. To make a conclusive 
statement on this possibility, one should study in detail the capability of LEP 2 to detect 
such a light Hi boson via its decays into r pairs, or into other hadronic modes. Most 
intriguingly, the set of parameters considered in this example also allow for a light right- 
handed stop squark and moderate mixing parameter, \X t \/Mq } of the type necessary to 
allow the possibility of electroweak baryogenesis |54[]. Hence, if such a Higgs boson cannot 
be discovered at LEP 2 via other decay channels, the final phase of LEP 2 will leave an open 
window for electroweak baryogenesis. A careful study of the CP-violating phases required 
for electroweak baryogenesis and the detection capabilities of LEP 2 for alternative decay 
modes becomes essential for testing this exciting scenario. 

We shall briefly comment on the enhanced LHC capabilities for Higgs-boson 



searches [pq| . The LHC has a considerably higher reach than LEP 2 and the upgraded 
Tevatron collider in the search for heavier Higgs bosons, and hence has more chances to 
unravel the complete Higgs-boson spectrum of the MSSM with explicit CP violation. At 
the LHC, Higgs bosons may be copiously produced via a wide variety of processes which 
depend in many different ways on the couplings of the neutral and charged Higgs bosons 
both to gauge bosons and fermions |5"5| . In the case of the CP- violating version of the 
MSSM under study, we have shown that mixing between states with different CP parities 
can dramatically modify those couplings and, hence, importantly affect the associated pro- 
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duction and decay mechanisms. Studies including CP-violating effects on the gluon-fusion 



production of Higgs bosons at the LHC have already been considered in the literature |56 
Our work provide the basic tool to improve further those studies, and to perform a com- 
plete analysis of the CP-violating effects on the many other Higgs-boson search mechanisms 
available at LHC. The LHC, together with the information gathered from experiments at 
LEP 2 and the upgraded Tevatron collider, will be capable of providing a thorough test of 
the MSSM Higgs sector and shed light on the possibility of explicit radiative breaking of 
CP invariance in supersymmetry. 

Finally, it is interesting to make a comparative analysis between our results and those 



obtained previously in [11]. In the latter work, the effective RG- improved potential was 
expanded up to renormalizable operators of dimension 4. The expansion was performed in 
powers of the stop-mass splitting, — m| , relative to the arithmetic average of the squared 
stop masses, M| USY = §(w~ + m \)- Moreover, the two- loop effect originating from the 
one-loop radiative corrections to the Yukawa couplings of the top and bottom quarks was 
not taken into account in the computation of the effective potential in [TT\\. Nevertheless, for 
moderate values of tan/3 and for all soft supersymmetry-breaking masses equal to M SUS y, 



the deviations of the results presented in JL1] with our results are expected to be small, 
for small values of the stop-mixing parameter X t = \A t — fi* cot /3|; the deviations will only 
grow for increasing values of X t . For larger values of tan/5, instead, the impact of the 
bottom-mass quantum corrections, which were omitted in pl| , is significant. In fact, only 
in the limit of small values of in which case the bottom-mass corrections are small, are 
both approaches guaranteed to give comparable numerical estimates. 

In Fig. |9|, we show the predictions for the mass of the lightest neutral Higgs boson 
H\ and its coupling to Z bosons, as obtained obtained by our RG approach (solid lines) 
and the operator-expansion method of |TT| (dashed lines). For the sake of comparison, 
we consider the same input parameters as those chosen in Fig. |], for vanishing gluino 
phase and three different values of the charged Higgs-boson mass: M#+ = 150, 200 and 
500 GeV. As was discussed above, we find that the predictions of the two works are in 
excellent agreement with one another for small values of X t . For large values of X t , instead, 
we observe larger quantitative differences in the results obtained by the two approaches, 
even though the qualitative behaviour of the two predictions exhibits a quite analogous 
functional dependence. 

Our one-loop RG-improved approach overcomes the limitations present in earlier 
analyses. In particular, our RG approach allows for a rather precise determination of 
the radiative effects on a generic Higgs-boson mass spectrum, even in cases of large stop 
mixings and/or large hierarchies between the left- and right-handed stop masses. Also, 
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within our RG approach, the important effects of the one-loop corrections to the quark 
Yukawa couplings are incorporated in the computation of the Higgs-boson masses and in 
their respective Higgs-boson couplings to gauge and fermion fields. A Fortran code that 
computes the Higgs-boson masses and couplings, including all the CP-violating effects as 
presented in this work, may be found in |57| . 



6 Conclusions 

We have performed a complete one-loop RG improvement of the effective Higgs potential 
in the MSSM, in which CP violation is induced radiatively by soft CP-violating trilinear 
interactions that involve the Higgs fields and the stop and sbottom squarks. Earlier studies 
P~0| , [TT] , [r2[ ] of the neutral Higgs-boson mass spectrum were based on a number of particular 
assumptions and/or kinematic approximations. The present work goes well beyond those 



studies, and extends the most detailed analysis [TT[, in which the one-loop RG-improved 
effective potential was expanded up to renormalizable operators of dimension 4, assuming 
a moderate mass splitting among the stop squarks. This assumption seemed to impose 
a serious limitation, given that CP-violating effects exhibit an enhanced behaviour for 
large values of the stop- mixing parameter X t . The results obtained using the present RG 
approach confirm, however, the qualitative phenomenological features found in ||11|| , for 
the Higgs-boson masses and their couplings to fermions and to the W and Z bosons. It 
offers very accurate predictions, at the same level as the most accurate calculations in the 



CP-conserving case (implying an uncertainty of order 3 GeV) p0|j43| , |44]| , for the Higgs- 
boson masses and for the whole range of the MSSM parameter space in the presence of 
non-trivial CP- violating phases. More specifically, the present study also includes two-loop 
leading logarithms associated with QCD effects and t- and 6-quark Yukawa couplings. It 
also contains all dominant two-loop non- logarithmic contributions to the one-loop effective 
potential, which are induced by one-loop threshold effects on the t- and 6-quark Yukawa 
couplings due to the decoupling of the third-generation squarks (as considered in [|43| . |44|j 
in the CP-conserving limit). These one-loop threshold terms, 5h u ^ and Ah UtC i, strongly 
depend on the phase of the gluino mass and so introduce new CP violation into the MSSM 
effective potential at the two-loop level. 

Large radiative effects of CP violation in the Higgs sector of the MSSM can have 
important phenomenological consequences on Higgs-boson searches at LEP 2, the Tevatron 
and the LHC. We have explicitly demonstrated that the radiatively-induced CP violation in 
the MSSM Higgs potential can lead to important effects of mass and coupling level crossing 
among the three neutral Higgs particles. These CP-violating effects of level crossing in 
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the Higgs sector modify drastically the Higgs-boson couplings to the up- and down-type 
quarks, and to the W and Z bosons. In particular, CP violation in the lightest Higgs 
sector becomes relevant for a relatively light charged Higgs boson, with Mjj+ < 160 GeV. 
For instance, for Mfj+ = 150 GeV and tan/? = 4, even a neutral Higgs boson as light as 
60 GeV may escape detection at LEP 2. However, the upgraded Tevatron may have the 
physics potential to explore such CP- violating scenarios at low tan/5 values, which may 
remain at the edge of accessibility even during the final LEP 2 run. 

We have also studied the effects induced by a non-trivial CP-odd gluino phase, which 
enters the effective potential at the two-loop level. The presence of a gluino phase gives rise 
to small but non-negligible changes in the Higgs-boson mass spectrum and in the couplings 
of the Higgs fields to the W and Z bosons. In this context, we have also found that the 
product of the scalar times the pseudoscalar coupling of the lightest Higgs boson Hi to 
the bottom quarks has a very strong dependence on the gluino phase. This product of 
couplings gives a measure of CP violation in the H\bb coupling, which can even be of order 
unity for relatively small charged Higgs-boson masses. 

Finally, it is worth stressing that CP violation decouples from the lightest Higgs sector 
in the large-mass limit of a heavy charged Higgs boson. This decoupling property, which 
was known to hold for the CP- violating self-energy effect [|T^,|rT|], has now been shown to be 
valid for the CP-violating vertex effects as well. As a result, the predictions for the lightest 
Higgs-boson mass and its couplings to gauge bosons in the above decoupling regime of the 
theory will be practically identical to the corresponding predictions in the CP-conserving 
case. However, unlike the lightest Higgs sector, CP violation does not decouple from the 
heaviest Higgs sector in the MSSM, opening up new possibilities for studying enhanced 
effects in CP-violating Higgs scalar-pseudoscalar transitions at the LHC or future muon 



colliders ]5T,52], where the heavy MSSM Higgs bosons can be resonantly produced. 

In conclusion: the present analysis has shown that the MSSM with explicit radiative 
breaking of CP invariance constitutes a very rich theoretical framework, introducing new 
challenges in the search for fundamental Higgs scalars at LEP 2, the Tevatron and the 
LHC. 
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Note Added 



After completion of the work described here, we saw which also computes the one-loop 
RG-improved effective potential of the MSSM with explicit CP violation. Here we further 
i) perform a more complete RG improvement of the effective potential, ii) develop a self- 
consistent treatment of the whole MSSM Higgs sector, taking into account the crucial one- 
loop relation between the charged Higgs-boson and neutral Higgs-boson mass matrices, and 
Hi) include the dominant two- loop non-logarithmic corrections to the effective potential, 
which may also have an important impact on the H \bb coupling. In the limits where a 



comparison between the results was possible, we find reasonable agreement with [58 
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A Derivatives of Background-Field-Dependent 

Masses 



In this Appendix, we list analytic expressions pertaining to derivatives of quark and squark 
masses with respect to their background Higgs fields. These derivative expressions are very 
useful, as they constitute the building blocks of the general one-loop effective potential 
presented in Section 2. We have divided the Appendix into three subsections. In the first 
subsection, we list the derivatives of quark masses with respect to Higgs fields, while the 
next two subsections contain the corresponding expressions for derivatives of squark masses 
with respect to neutral and charged Higgs fields, respectively. 



A.l Quark Derivatives 

First, we give the derivatives related to non-vanishing tadpole contributions: 

1 /*»?\ - |M>, = (A.1) 



V 2 \ OCp2 I V\ \ 0<pi 

where the operation ( • ■ • ) denotes that the above expressions should be evaluated in the 
ground state of the Higgs potential. 

Then, the self-energy-type derivatives involving neutral and charged Higgs bosons 
may be listed as follows: 

5 2 m?\ /d 2 m 2 \ ,, l9 /d 2 m 2 \ I d 2 m 2 



d4>l I \ da 2 . I \ 50 2 I \ da 2 

I d 2 m 2 \ \h b \ 2 m 2 / d 2 m\ \ \h b \ 2 ml 



\ d(j)i d(p 1 I m 2 — ml \ d(fil d(j) l I m 2 h — m 2 1 



d 2 m\ \ \h t \ 2 m 2 I d 2 ml \ \h t \ 2 m 2 



j 2 u V2 i m 2 t -m 2 b \ d(pt d(p 2 I m 2 b - m 2 t ' 

d 2 mf \ I d 2 m\ \ I d 2 m 2 h \ I d 2 m 2 h \ \h t h b \m t m b 



(A.2) 



, d(pi (902 / \ d(pt 50x / \ 50i" 502 / \ 502~ 50x / ml — m 2 

Note that we have not listed derivatives that vanish. 
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A. 2 Derivatives of Squark Masses with Respect to 
Neutral Higgs Fields 

This section contains the derivatives of the field-dependent squark masses m 2 and m 2 ^ . 
with respect to neutral Higgs fields 4>i,2 an d 01,2- We first give the tadpole terms: 



1 / dm 



ti(t 2 ) 



vi \ d(f)i 



1 / dm l(b 2 ) 



vi \ d(pi 



1 /dm 2 



ti(t 2 ) 



v 2 \ dfa 



1 / dm l(b 2 ) 



V 2 \ d(p2 



1 / dm 2 



h(t 2 ) 



v 2 \ dai 
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v 2 \ dai 



gl + g' 2 
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x t 
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x 6 (Mq - M fe 2 - ±x v 2 cos 2 
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2 1 
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(A.3) 



where the coupling parameters x t and x are defined after ( j2.13j ). Then, the self-energy-type 
terms (d 2 m 2 k /daidaj), with q k = t 1 , 6 1( £ 2 , &2 an d z,j = 1, 2, are found to be 



d2m i(t 2 ) 
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, d2 K(t 2 ) 

i dai da 2 

d2m l(b 2 ) 
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(A.4) 



6i 62 ' 

In addition, the non- vanishing CP- violating self-energy terms (<9 2 m; 2 / 'dfadcij) are given by 
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(A.5) 



Finally, the CP-conserving self-energy- type derivatives (<9 2 m 2 fc /<90j<90j) have been calcu- 
lated to be 
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A. 3 Derivatives of Squark Masses with Respect to 
Charged Higgs Fields 

Here we evaluate the derivatives of the field-dependent squark masses with respect to 
charged Higgs fields. To calculate the expressions (d 2 m 2 k / dtfif d(f>~) directly turns out to 
be a formidable task. The reason is that m 2 k are the eigenvalues of a non-trivial (4 x 4) 
squark mass matrix M 2 (cf. (|2.11 )), and their analytic form is very complicated. Therefore, 
we proceed differently, using a mathematical trick which was first applied in [pSfl . 

First, we notice that (dm 2 k /d(j)f) = 0, as a consequence of the fact that the true 
ground state of the effective potential should conserve charge. Then, one may make use of 
the eigenvalue equation: 



det (M< 



m l + Am l + Bm l + Cm i + D 







(A.7) 



with 
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c 

D 



~ (TrlM 2 - TrM & 



to obtain 
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A a m l k + B a m \ k + + D u 



n o 2 



(A.9) 
etc.. For our 



with qi,q k = ti,bi,t 2 ,b 2 , and A {j = d 2 A/d(j)id(j)j, B ii = d 2 B /d(f)Jd(j)J , 
purposes, it is sufficient to calculate the derivatives with respect to (fif and 4> 2 . m particular, 
it proves convenient to use a representation in which the columns and rows '2' and '3' of 
the (4 x 4) matrix M 2 have been interchanged. With such a reordering, we find 
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where M\ and Ml are the usual t- and 6- (2 x 2)-mass matrices, respectively, and 
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Then, the relevant coefficients A 12 , -B12, C12 an d -D12 may be expressed in a compact form 
as follows: 
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(A.12) 



With the help of ( |A.12j ), it is straightforward to obtain the derivatives d 2 m 2 k /d<f)fd<j) 2 ■ 
More explicitly, we have 
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B Higgs-Boson Masses and Mixing Angles 



Here we present analytic expressions for the Higgs-boson masses M H .(m t ) (i = 1,2,3) 
and the corresponding (3 x 3) orthogonal matrix O, after diagonalizing the RG-improved 
Higgs-boson mass matrix (M 2 N )(m t ). 

For notational simplicity, we do not display explicitly the functional dependence of 
M% on m t . The mass eigenvalues of the (3 x 3) matrix M% are then obtained by solving 
the characteristic equation of cubic order: 

x 3 + rx 2 + sx + t = 0, (B.l) 

with 

r = -Tr(M 2 N ), 
1 



S = 2 



Tr 2 (M 2 N ) - Tr(M%) 
t = -det(M%). (B.2) 



To this end, it proves useful to define the following auxiliary parameters: 

3s - r 2 

P = 



3 

2r 3 rs 

9 = 1 ~ 7 + t ' 

= P_ + Q_ (B.3) 
27 4 K ' 

To ensure that the three eigenvalues are positive, it is necessary and sufficient to require 
that 

D < 0, r<0, s > 0, t < 0. (B.4) 

Imposing these inequalities on the kinematic parameters of the theory, we may express the 
three mass eigenvalues of M 2 N as 



= - l -T + 2v /^/3cOs(| + y), 

= - l -r + 2v /^/3cos(| - y), 

= -\r + 2 v / I p/3cos(|), (B.5) 
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with 



arccos 



Q 



and 



< If < 7T. 



;b.6) 



v 2^-^/27' 

Since the Higgs-boson mass matrix Ai% is symmetric, we can diagonalize it by means of an 



orthogonal rotation O as stated in (|3.30|) . Furthermore, one can show [59] that the Higgs- 
boson mass eigenvalues in (|B.5|) satisfy the desired mass hierarchy in accordance with the 



inequality of (3.31 



If M?j, with i, j = 1,2,3, denote the matrix elements of the elements Oij can then 
be obtained by appropriately solving the underdetermined coupled system of equations, 
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More explicitly, we have 
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(B.10) 



In ( |B.10[ ) , the abbreviation s x = sign (x) is an operation that simply gives the sign of a 
real expression x. 
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Figure 1: Numerical estimates of (a) Mu x and Mh 2 and (b) g\ iZZ as functions o/arg (A t ), 
for the indicated choices of MSSM parameters. Solid lines correspond to arg (m,g) = 0, 
dashed lines to arg (m g ) = 90°. 
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Figure 2: As Fig. g, but with n = 4 TeV, \A t \ = \A b \ = 2 TeV and M SU sy = 1 TeV. 
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Figure 3: As in Fig. ^ but with tan/5 = 20. 
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Figure 4: As in Fig. |, but with M H + = 135 GeV. 
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Figure 5: Numerical estimates of (a) M Hl and (b) 2\M H2 - M H3 \/{M H2 + M H:i ) as func- 
tions of the CP-violating phase arg(A t ). Lower values of the same line type correspond to 
arg (m,g) = 0, the higher ones to arg (rrig) = 90°. 
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Figure 6: Numerical estimates of (a) (gf Ilbb ) 2 + (g^bb) 2 and (b) 1\(g s Hibb ) | /[(^ lb6 ) 2 + 

(g^bb) 2 } versus arg(A t ). Solid lines correspond to arg (rrig) = 0, dashed ones to arg (m,g) = 
90°. 
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Figure 7: As in Fig. [6|, but with tan/5 = 20. 
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arg (A) = arg (A h ) [ deg ] 

(a) 




Figure 8: Numerical estimates of (a) Mu x and Mh 2 and (b) [{gfj^) 2 + (fi^M,) 2 ] versus 
aig(A t ), with soft squark masses: Mq = M b = 0.6 TeV and M t = 0. Solid lines are for 
arg (rrig) = —90° and dashed lines for arg (rrig) = 0°. 
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Figure 9: Numerical predictions for Mh 1 and gjj^z as obtained by the present complete 
RG approach (solid lines) and the operator- expansion method of / [Ji|/ (dashed lines). 
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